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ABSTRACT 

Type la supernovae (SNe la) are thought to originate from the thermonuclear explosions of 
carbon-oxygen (C-O) white dwarfs (WDs). The gravitationally-confined detonation (GCD) model 
is a well-explored explosion model where unstable thermonuclear burning initiates in an accreting, 
Chandrasekhar-mass WD and forms an advancing flame. Rising, burning material breaks out of the 
WD surface, engulfs the WD, and leads to a gradient-induced detonation in the collision region, which 
subsequently unbinds the WD. We show that when the flame consumes enough fuel, the conditions 
for detonation are not met, and the WD does not detonate. We show through 2D and 3D simulations, 
such failed-detonation SNe expel a few 0.1 M of burned and partially-burned material, but do not 
unbind the WD. A fraction of the material falls back onto the WD, polluting the remnant WD with 
intermediate-mass and iron-group elements, that likely segregate to the core forming an iron-core C-0 
WD. The remaining material is asymmetrically ejected at velocities comparable to the escape velocity 
from the WD, and in response, the WD is kicked to velocities of a few hundred km s^ 1 . Such kicks 
may unbind the binary and eject a runaway/hyper- velocity WD. Although the energy and ejected 
mass of the failed-detonation SN are a fraction of typical thermonuclear SNe, they are likely to ap- 
pear as sub-luminous low-velocity SNe la. Such failed-GCDs might therefore explain or are related 
to the observed branch of peculiar SNe la, such as the family of low-velocity sub-luminous SNe (SN 
2002cx/SN 2008ha-like SNe). 

Subject headings: supernovae: general — supernovae: individual (2002cx, 2008ha) — hydrodynamics 
— white dwarfs — ISM: supernova remnants 



1. INTRODUCTION 

Type la supernovae (SNe la) are among the most ener- 
getic explosions in the known universe, releasing ~ 10 51 
ergs of kinetic energy in their ejecta, and synthesiz- 
ing ~ 0.7 Mq of ra dioactive 56 Ni . The discovery of 
the Phillips relation (jPhillipsI 119930 enabled the use of 
SNe la as standardizable cosmological candles, and has 
ushered in a new era of astronomy lead ing to the discov- 
ery of the acceleration of the universe (jRiess et al.lll998t 
iPerlmutter et alJll999t l. leading to the 2011 Nobel Prize 
in physics. 

In the following we introduce an unexplored variant of 
well-studied single-degenerate (SD) model of SNe la in 
which a Chandrasekhar-mass, carbon-oxygen (C-O) WD 
accretes mass from a non-degenerate main sequence or 
red giant companion. As the WD approaches the Chan- 
drasekhar limit, carbon burning is initiated in its con- 
vective core. After a few hundred years of this "simmer- 
ing" phase, unstable thermonuclear burning is expected 
to ignite at one or more off-centered points, generating 
a buoyantly-rising, subsonic deflagration flam e bubble 
(|Kuhlen et alJl200a IZhTgale et al.ll200l I20T1 . We re- 
fer to this phase as the "deflagration phase" . 
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We adopt the assumption made by the gravitationally- 
confined detonation (GCD) scenar i o (Dordan et alj|2008t 
iMeakin et aTJ 120091 : Uordan et all |201_2|) of SNe la in 
which no transition to detonation occurs as the ris- 
ing flame approaches the low-density surface layers of 
the WD as is posited by the "deflagrati on-to-detonation 
transition" (DDT) model of SNe la (|Khokhlovl U99H; 
IGamezo et all l200l 120051) . We assume an off-centered 
distribution of ignition points as in the GCD scenario. 
We further focus on the "puls ationally-assisted" variant 
of the GCD model described in lJordan etafl (pfMIMl 
which we review in section [2] 

We demonstrate that under these assumptions, if 
enough mass is burned during the deflagration phase, the 
conditions to trigger a detonation via the GCD mech- 
anism are not realized and the WD fails to detonate. 
Such "failed-detonation" (FGCD) models have numer- 
ous remarkable implications for the observable proper- 
ties of the resulting explosion and its outcomes. These 
include the production of peculiar SNe la events with 
low expansion velocities, low luminosities and low ejecta- 
mass; all broadly consistent with the observed proper- 
ties of a branch of peculiar SNe Iasimilar to SN 2002cx 
and/or SN 2008ha. Even more remarkably, we demon- 
strate that the WD survives the initial la event, receiv- 
ing a large velocity kick from the asymmetric nature of 
the initial deflagration, and is significantly enriched with 
both intermediate-mass (IME's) and iron-group elements 
(IGE's), forming a peculiar heavy/iron core WD. 

2. PHYSICS OF THE FAILED-DETONATION SCENARIO 

The pulsat i onally -assisted GCD model described in 
I Jordan et ail (|2012t ) be gins with an off-centered ignition 
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region in a Chandrasekhar-mass WD which forms an 
advancing thermonuclear flame. The hot, buoyant ash 
quickly rises forming a large plume that advances and 
burns towards the surface of the star. Energy imparted 
to the WD during the deflagration causes the star to 
pulse. As the WD expands, hot ash erupts from the star, 
flows around the stellar surface, and mixes with cold fuel. 
The mixture collides at the antipodal point on the sur- 
face from breakout. Meanwhile, the WD reaches max- 
imum expansion and contracts. As the star contracts, 
it squeezes the fuel-ash mixture in the collision region 
to temperatures and densities that considerably shorten 
the minimum length scale on which the Zel'dovich gra- 
dient mechanism can trigger a detonation. A detonation 
is triggered in the mixture which subsequently consumes 
the WD. 

Roughly speaking, typical fuel-ash mixture properties 
associated with an induction gradient < 10 km are den- 
sities >lx 10 7 g /cm 3 and peak temperatures > 1 x 10 9 
K (seelArnett fc Livndll994HNiemever fe Wooslevlll997t 
Khok hlov et al.ll997tlRopke et alj2007HSeitenzahl et al.l 
20093, for a complete discussion). 

The FGCD scenario proceeds similarly to the pulsa- 
tional GCD model; however, the scenarios deviate in the 
amount of burning — and therefore the amount of energy 
released — during the deflagration phase. More burning 
occurs during the deflagration phase of the FGCD than 
in the GCD which means more mass is consumed by 
the flame and ejected from the star in the form of ash, 
and more energy is delivered to the WD. The WD is 
thus modified to a higher degree than in the GCD model 
and as a result, the WD does not attain high densities 
and temperatures in the collision region during its con- 
traction. The fuel-ash mixture never reaches the critical 
conditions for detonation and therefore the star never 
detonates. 

The deflagration is strong enough to prevent the for- 
mation of the conditions for detonation, but too weak to 
completely unbind the star. The WD will thus remain 
intact; though, it will have a lower mass. Interestingly, 
since the ignition of the deflagration is off-centered, the 
remaining bound stellar material is kicked by the ejec- 
tion of the ash and obtains a velocity of hundreds of km 

This system produces an asymmetric outburst of a few 
tenths of a solar mass of deflagration products rich in 
IME's (such as Mg, Si, and S) and iron-group elements 
IGE's (such as Fe, Co, and Ni). Since the deflagration 
only processes a few tenths of a solar mass of material, 
the FGCD produces a small total mass of 56 Ni compared 
to normal SNe la which convert nearly half of the WD 
to 56 Ni. 

Some of this material attains escape velocity and some 
falls back onto the star. The velocity of the outflow will 
be slow (approximately a few thousand km s _1 ) related 
to normal SNe la due to the comparatively small amount 
of energy released in the FGCD scenario. 

3. SIMULATIONS OF THE FGCD MODEL 

3.1. Simulation Setup 

We used the Adaptive Mesh Refinement ( AMR) 
FLASH applicati on framework (|Dubev et al.l l200a 
iFrvxell et al.l l2C)00) to perform our simulations of FGCD 



models. FLASH has been previously used to simulate 
SNe la in both 2D cylindr i cal and 3D Cartesia n geometry 
(jJordan et al.ll20H [20081 iMeakin et al.ll2009ft . Our sim- 
ulations include an advection-diffusion- rea ction (ADR) 
treatment of the thermonuclear flame dCalder et al.1 
[2007t iTownslev et al.l [20071 : iSeitenzahl et al.ll2009bh . an 
equation of state that includes contributions from black- 
body radiation, ions, and electrons of an arbitrary degree 
of degeneracy (T immes fe Swestvll2000l ). and the multi- 
pole t reatment of gravity includ ed in the FLASH stan- 
dard (jASC FLASH Centeril20"Il . 

We performed four 3-dimensional (3D) exploratory 
simulations to test the feasibility of the FGCD and 1 
full 2-dimensional (2D) model to adequately observe the 
FGCD at late times. 

Our 3D simulations included a reduced domain size 
and a moderate resolution (8 km) which reduced their 
computational expen se. We initiali z ed the se simulations 
similarly to those in Uordan et al.l (|2012t ) with a 1.365 
Mq WD placed at the origin of the domain. We used 
the same distribution of sixty three, 16 km radius igni- 
tion "points" distributed in a 128 km radius spherical 
volume. We chose 48 km, 38, km, 28 km, and 18 km as 
offset distances along the z-axis of the spherical volume. 
Table[T]lists the initial conditions for each simulation and 
the corresponding names we gave them. We ran the sim- 
ulations from ignition, through peak stellar expansion, 
and at least until the WD reached a maximum central 
density upon contraction. 

Our 2D simulation was performed with a large domain 
in 2D cylindrical geometry at 4 km resolution and ran for 
60 seconds. This simulation was initialized similarly to 
the 3D simulations except we placed only 4 bubbles in a 
64 km spherical volume offset by 70 km along the z-axis 
(axis of symmetry). We chose these initial conditions to 
obtain an FGCD in 2D given what we learned from our 
3D simulations. The larger domain size allowed us to 
follow the outer layers of the ejecta for the entirety of 
the simulation. 

3.2. Simulation Results 

We obtained an FGCD from each of our simulations, 
and were able to quantify features discussed in section O 
Table [1] contains a collection of these quantities from our 
simulations. 

The energy liberated during the deflagration is between 
89% and 150% of the binding energy of the WD. How- 
ever, in each simulation the deflagration fails to unbind 
the WD. The WD expands in response to the deflagra- 
tion, reaches its maximum level of expansion, and then 
contracts. Figure Q] shows the evolution of the central 
density, p c , of the WD and illustrates the oscillatory na- 
ture of the WD after the deflagration. The more en- 
ergy released during the deflagration, the more the star 
expanded, the longer the pulsational period, and the 
smaller the p c at maximum contraction. 

The fact that the deflagration releases more than 100% 
WD binding energy in some models suggests that the 
WD should be unbound after the deflagration; however, 
the entire energy budget does not work to only unbind 
the star. For example, some of the energy is lost when 
the ejecta and the remnant WD are accelerated to high 
velocities. Thus, even though the deflagration liberates 
enough energy to completely unbind the star, the energy 
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TABLE 1 

Simulations Properties 



sim Aa; a n ign b r ign c z ign d E nuc <= C-0 B f IME B IGE B v B s C-0 E h IME E IGE E E^.e 1 v avc ,! 

name (km) (km) (km) (E Mn ) (M ) (M Q ) (M ) (km s' 1 ) (M ) (M ) (M ) (ergs -10 50 ) (km s" 



2D70 4 4 64.0 70.0 0.89 0.93 0.07 0.13 119 0.13 0.03 0.07 0.32 3,731 

3D48 8 63 128.0 48.0 1.06 0.84 0.06 0.09 351 0.16 0.05 0.16 0.90 4,93: 

3D38 8 63 128.0 38.0 1.21 0.71 0.07 0.08 383 0.24 0.06 0.20 1.3 5,09! 

3D28 8 63 128.0 28.0 1.39 0.56 0.05 0.09 520 0.32 0.09 0.25 1.8 5,22: 

3D18 8 63 128.0 18.0 1.53 0.43 0.04 0.08 406 0.43 0.10 0.28 2.2 5,21: 



a Maximum spatial resolution. 

Number of ignition points. 
c Radius of the spherical volume containing the ignition points. 

Z coordinate of the origin of the spherical volume containing the ignition points. 
e Energy released during the deflagration phase divided by the binding energy of the star (E^in — 4.5 X 10 50 ergs). 

"B" refers to gravitationally bound material. 
g Velocity of the gravitationally bound material . 

"E" refers to material that will escape the system. 
1 Average kinetic energy of the escaping material. 
J Mass weighted velocity of escaping material. 
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Fig. 1. — The central density, p c , of the WD vs time. Note that 
the p c achieved during contraction is < 10 7 g/cm 3 , the density at 
which the Zel'dovich gradient mechanism can produce a detonation 
on length scales less than ~ 10 km. 

is partitioned in such a way that a portion of the original 
WD remains gravitationally bound. 

In each of our models, the WD gets a kick in the op- 
posite direction from which the buoyant ash rises and 
breaks through the surface of the star. We measured this 
velocity to be on the order of hundreds of km s _1 and list 
these values for each simulation in table [T] 

Some of the material from the FGCD escapes and 
achieves high velocities, while some of the material is 
bound to the remnant WD and will eventually accrete 
onto its surface. Table Q] lists the composition of mate- 
rial that remains gravitationally bound and eventually 
will mix with the remnant WD. In all cases between 0.1 
M Q and 0.2 M of IME's and IGE's (some of which is 
radioactive 56 Ni) remained bound to the star. The evolu- 
tion of the WD as the material falls back onto and heats 
the star is an interesting question and one which we will 
examine in future work. 

The composition of the material that escapes the sys- 
tem is also listed in table Q] along with the kinetic energy 
and the mass-weighted velocity of the ejecta. This mate- 
rial includes carbon and oxygen, IME's, and IGE's, and 



ranges from 0.2 to 0.8 M®. In general, the more material 
that is burned during the deflagration, the more material 
that escapes. 

Figure [2] shows shows the density structure and com- 
position with overlaid velocity contours of the 2D model 
at 60s. Note that the density profile of the FGCD model 
is asymmetric in velocity space between the hemisphere 
corresponding to the ejected deflagration and the oppo- 
site hemisphere of the system. The remnant of the WD 
can be seen as the tiny high-density feature slightly below 
the origin of the domain. 

The figure also shows the nature of the asymmetry in 
composition of the structure. The north side of the rem- 
nant contains the products of the deflagration that were 
sprayed from the surface of the WD. Clumpy structures 
of IME's and IGE's exist at a range of velocities in the 
northern hemisphere of the domain but are less abundant 
in the southern hemisphere. The asymmetries suggest 
that this object would look much different depending on 
the viewing angle of the observer. 

All of our models produce a relatively small amount 
of radioactive 56 Ni. Though we do not perform detailed 
nucleosynthetic post processing of the 3D models, we can 
set an upper limit on the 56 Ni yields with the amount of 
IGE's produced, which ranges from 0.2 M Q to 0.36 M . 
Neutronization through electron capture reactions dur- 
ing the deflagration would shift production away from 
56 Ni though and reduce its contribution to the IGE to- 
tals. Whether even lower masses of IGE's (and thus 
56 Ni) could be produced under appropriate conditions 
(e.g. comparable to that observed in SN 2008ha, a very 
faint SN with extremely low ejecta velocities), is yet to 
be explored. 

4. DISCUSSION AND PREDICTIONS 

4.1. Sub-luminous Low-Velocity SNe 

The most prominent features of the FGCD models are 
their low-mass, low-velocity ejecta and overall low ki- 
netic and nuclear energy. These properties translate into 
the production of typically sub-luminous, low-velocity 
SNe la. It is therefore natural to examine whether SNe 
with such characteristics have already been discovered. 
In particular, one may explore peculiar SNe exhibiting 
either extremely low mass ejecta, such as SN 2008ha 
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Fig. 2. — Images from the 2D70 simulation. The black contour marks the transition between gravitationally bound and escaping material. 
The white contours are of the magnitude of the velocity field. From the inner most contour and moving outwards, the values are: 2,000 km 
s _1 , 4,000 km s _1 , 6,000 km s" 1 , 8,000 km s" 1 , and 10,000 km a -1 . The images show the simulation at 60 seconds, (left) Log density of 
the remnant. The values of the density are given by the color bar on the left, (right) The composition of the remnant. Blue is C-O, green 
is IME's, and red is IGE's. 



(Foley et all I2009t IValenti et al 
SNe ([Perots et all iSoToT 



[20(51 SN 2 005E-like 
Kasliwal et all I2012al). or SN 



2002bi-like SNe (iKasliwal et al.ll2012bt [Poznanski et al.l 
l2Tjl0t berets et al.ll2011l). or low e j ecta velocity s uch a s 
SN 2002cx-like SNe lLi et al.l (l2Tj0l) : lBrarTcli et al.l d200l : 
Uha et all (|2006l ). Some of the for mer low-mass-ejecta 
SNe, however, may present He signature in their spectra 
(type lb SNe), inconsistent with our models. We there- 
fore focus on observed low velocity, faint SNe la. 

Normal type la SNe differ in their ejecta velocities as 
measured in some standard method, but they generally 
fall between 9000 — 14000 km s _1 near peak luminosity, 
with simi lar dispersion at lat er times (as derived from the 
Si II line lBender et al.ll2005 | ). The velocit i es of even the 
lowest velocity SNe in the iBender et ail (|2005f ) sample 
much exceed the mass-averaged FGCD velocities listed 
in table [TJ We also note a trend of more energetic and 
likely more luminous (larger IGE yield) SNe to be ac- 
companied by higher ejecta-velocities over almost an or- 
der of magnitude in kinetic energy. The only other type 
of SNe with such low expansion velocities are the branch 
of peculiar type la SNe, na med for the pro t otype for this 
class of supernovae, 2002cx ([Li et al.ll2003t iBranch et al.l 



[20TIl lJha~e"t alJ2006|): such SNe may a l so hav e an enregy- 
velocitv correlation IMcClelland et ahl ([2010D as observed 
in our simulations. 

SN 2002cx-like events are characterized by luminosities 
which lie too low in comparis on to the Phill ips relation 
for Branch-normal l a events dLi et al.l [2003) , low pho- 
tosphcric velocities ([Li et al.l 120031 ). weak intermediate- 
mas s element lines (inc luding Si II, S II, Si III, and Ca 
II) ([Branch et al.l [2004^ . and late-time o ptical nebular 
spectra dominated b y narrow Fe II lines (IBranch et al.l 
[200llJha et^[20TJl . Since the discovery of SN 2002cx, 
a number of other 2002cx-like events h ave been discov- 
ered including 2002es, 2005P, 20 05hk ([Chornock et al.1 
MM), 2008ge (F oley et al.l 120 lOal ) and 2008ha. The lat- 
ter event (SN 2008ha), in particular, is consistent with 
extremely low mass ejecta and energetics. We predict 
the FGCD models to produce similar properties to those 
characterizing SN 2002cx like SNe, given the low expan- 
sion velocities and the low estimated 56 Ni yield, and po- 
tentially even explaining SN 2008ha like events with low 
mass ejecta. 

Though our initial set of simulations is limited, the ro- 
bust features of FGCD's, including low velocity ejecta, 
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the expected low luminosity (due to the small yield of 
56 Ni) and their low mass ejecta (comparable to that SN 
2008ha) make them tantalizing candidate progenitors for 
this branch of peculiar SNe. Note that the single degen- 
erate origin of such SNe is also consistent wit h the overall 
typically young (but not necessarily young; iFolev et al.l 
2010) environments found for SN-2002cx like SNe, com- 
pared to the expectations from, e.g. double degenerate 
SNe (older environments) or core-collapse SNe (only very 
young environments). 

4.2. Heavy metal/Iron- core WD 

In our FGCD scenario, a large amount of burnt mate- 
rial falls back to the remnant WD. From table [T] the WD 
may incorporate as much as 0.1 Mq of IGE's and 0.07 
Mq of IME's of fallback material, together comprising 
as much as ~ 1/4 of the remnant C-0 WD. Such heavy 
elements are likely to gravitationally settle, in time, to 
the WD core, making these heavy-metal polluted WD's 
become iron/heavy-core WDs. The existence of iron- 
core WD's has been considered be fore, with even the po- 
tential observation of such WDs (jProvencal et al.l Il998b 
iCatalan et aHl2008l ). The FGCD scenario therefore pro- 
vides a novel evolutionary scenario for the formation of 
these iron/heavy-core C-0 WD's. A somewhat related 
scenario of failed SN was sugge sted for the forma tion of 
O-Ne-Mg WDs with iron cores (|Isern et al.lll991l ). 

4.3. WD Natal Kicks 

FGCD's produce a highly asymmetric ejection of ma- 
terial. This is not unique amongst various models for 
SNe explosions. However, in our FGCD case, the WD 
survives the explosion. Considering momentum conser- 
vation, this gives rise to a unique outcome, namely that 
the surviving WD is kicked at very high velocities, rang- 
ing hundreds of km s _1 . The FGCD scenario suggests 
the existence of strong WD natal kicks, and provides 
an interesting prediction per the existence of hyperve- 
locity WDs. Taken together, the potential existence of 
a a heavy core WD (discussed in section 14.21) , and the 
high ejection velocity produce a highly peculiar object, 
which, if observed may provide a possibly unique smok- 
ing gun signature. One should note, however, that the 
population of halo WDs may also have relatively high 
velocities, and it might therefore be difficult to pinpoint 
the kinematic property as related to a natal kick (unless 
the WD is massive and young; an unlikely possibility for 
the old population of halo WDs) . 

We note that velocities of hundreds of km s~ x are much 
larger than the orbital velocities of most binaries, and a 



kick velocity of such magnitude is likely to unbind the 
binary (ejecting the companion star at high velocity, but 
much lower than the WD kick velocity). WD's accreting 
from evolved stars could have a large binary separation, 
and be disrupted following the WD kick. However, WD's 
accreting from a main sequence companion might have 
highly compact configurations. In this latter case, the 
WD kick might not unbind the binary, and may only 
change its orbit, and provide an overall kick to the binary 
center of mass as a whole. An interesting aspect is then 
the result of a following accretion epoch involving the 
iron-core WD, and its possible evolutionary outcomes. 
Finally, in both cases of evolved and main sequence mass 
donors, the WD natal kick may eject the remnant WD 
towards the companion, in which case a collision may 
occur, with potentially interesting explosive results. 

5. SUMMARY 

We have discovered a scenario in which the standard 
GCD model fails to detonate due to a large amount of 
burning during the deflagration phase. The failed GCD 
results in an asymmetric outburst of deflagration mate- 
rial consisting of intermediate elements and iron group el- 
ements along with a fraction of the original WD still grav- 
itationally bound. The models produce a faint SNe la 
with a slowly evolving light curve due to the low 56 Ni 
yield and the low energetics. The remaining WD gets a 
kick on the order of hundreds of km s _1 and is contami- 
nated with fall-back from the deflagration, producing an 
iron/heavy core WD. We presented our initial simula- 
tions to quantify the some of the bulk observable proper- 
ties and hypothesize that the FGCD model is a possible 
explanation for 2002cx-like SN. Future studies will ex- 
plore the detailed observational features of FGCD SNe 
and their direct comparison to observations. 
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